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Frequency locking by external forcing in systems with rotational symmetry 

Lutz Recke^^, Anatoly Samoilenko-f ^, Viktor Tkachenko^^, and Serhiy Yanchuk^"!! 

Abstract. We study locking of the modulation frequency of a relative periodic orbit in a general S^-equivariant 
system of ordinary differential equations under an external forcing of modulated wave type. Our 
main result describes the shape of the locking region in the three-dimensional space of the forc- 
ing parameters: intensity, wave frequency, and modulation frequency. The difference of the wave 
frequencies of the relative periodic orbit and the forcing is assumed to be large and differences of 
modulation frequencies to be small. The intensity of the forcing is small in the generic case and can 
■ be large in the degenerate case, when the first order averaging vanishes. Applications are external 

electrical and/or optical forcing of selfpulsating states of lasers. 
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1. Introduction. This paper concerns systems of ordinary differential equations of the 
C/j dx 

Q: — = fix)+^g{x,/3t,at). (1.1) 

^ ; Here / : M'' and g : x R*^ are C'-smooth with Z > 3, and a > 0, /3 > 0, 7 > 

■ are parameters. The vector field / is supposed to be 5^-equivariant 

fie^^x) = e^^fix) (1.2) 

I for all X G R" and ^ € R, where ^4 is a non-zero real n x n-matrix such that = —A and 

(<~^ ■ e'^'^^ = I. The forcing term g is supposed to be 27r-periodic with respect to the second and 

Q> ■ third arguments 

^ ! gix, tp + 27r, ip) = g{x, V', (f) = g{x, V', + 27r), 

QQ ■ and to possess the following symmetry 

^" g{e^^x,i;,^ + i) = e^^g{x,i;,^) (1.3) 

for all X G R" and ip,ip,(, G R. Finally we assume that the unperturbed system 

K. (1.4) 
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has an exponentially orbitally stable quasiperiodic solution of the type 

x{t) = e^"«*xo(/?oi). (1.5) 

Here : K — )■ M" is 27r-periodic, and > 0, /3o > are the wave and modulation frequencies 
of the solution (1.5). The main property of quasiperiodic solutions of the type (1.5) is that 
they are motions on a 2-torus without frequency locking. Those solutions are sometimes 
called modulated waves [16] or modulated rotating waves [5] or relative periodic orbits [9]. 
We assume that the following nondegeneracy condition holds: 

The vectors Axq{'iP) and "j^i'^P) linearly independent. (1-6) 

It is easy to verify that (1.6) is true for all -i/^ G M if it holds for one ijj. Condition (1.6) implies 
that the set 

:= {{e^^xoi^)) G : v', V G K}, (1.7) 

which is invariant with respect to the flow corresponding to (1.4), is a two-dimensional torus. 

Our main results describe open sets in the three-dimensional space of the control param- 
eters a, (3 and 7 with |a — ao| S> 1 and (5 k, I5q, where stable locking of the modulation 
frequencies /3 of the forcing and /3o of the solution (1.5) occurs, i.e. where the following holds: 
For almost any solution x[t) to (1.1), which is at a certain moment close to Ti, there exists 
(7 G M such that 

inf \\x{t) - e'^'^xo{l3t + a)\\ ^ for large t. (1.8) 

Our results essentially differ in the so-called non-degenerate and degenerate cases (see 
Section 2). The degenerate case includes all cases when the averaged (with respect to the fast 
wave oscillations) forcing 

1 /■^'^ 

fi(i(x,V') := — g{x,^,ip)dip (1.9) 

vanishes identically, for example, if g is of the type g{x,ip,ip) = e"^'^^('0). Roughly speaking, 
under additional generic conditions the following is true: 

Non-degenerate case: If a is sufficiently large and 7 is sufhciently small and j3 — (3q is of 
order 7, then locking occurs. 

Degenerate case: If a is sufficiently large and 7^/a is sufficiently small and 13 — (Sq is of 
order 7^/a, then locking occurs. 

The present paper extends previous work on this topic for particular types of the vector 
field / and the forces g [23, 20, 19]. In particular, the considered special cases in [23, 20, 19] 
are doubly degenerate in the sense that not only the averaged forcing (1.9) vanishes but also 
the second averaging term turns to zero. Remark that in [18] related results are described for 
the case that both differences of modulation and wave frequencies are small, and [17] concerns 
the case when the internal state as well as the external forcing are not modulated. For an 
even more abstract setting of these results see [4]. 

Systems of the type (1.1) appear as models for the dynamical behavior of self-pulsating 
lasers under the influence of external periodically modulated optical and/or electrical signals. 
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see, e.g. [15, 1, 12, 13, 14, 24, 26], and for related experimental results see [8, 22]. In 
(1.1) the state variable x describes the electron density and the optical field of the laser. In 
particular, the Euclidian norm describes the sum of the electron density and the intensity 
of the optical field. The S'^-equivariance of (1.4) is the result of the invariance of autonomous 
optical models with respect to shifts of optical phases. The solution (1.5) describes a so-called 
self-pulsating state of the laser in the case that the laser is driven by electric currents which 
are constant in time. In those states the electron density and the intensity of the optical field 
are time periodic with the same frequency. Self-pulsating states usually appear as a result of 
Hopf bifurcations from so-called continuous wave states, where the electron density and the 
intensity of the optical field are constant in time. 

External forces of the type ^e^"^g{f3t) appear for describing an external optical injection 
with optical frequency a and modulation frequency /3. In spatially extended laser models 
those forces appear as inhomogeneities in the boundary conditions. After homogenization of 
those boundary conditions and finite dimensional mode approximations (or Galerkin schemes) 
one ends up with systems of type (1.1) with general forces of the type "fgix, /3t, at) with (1.3). 
External forces of the type 

g{x,/3t, at) = gi{x, /3t) with gi (e^^x, tp) = e^^gi (x, ip) 

appear for describing an external electrical injection with modulation frequency (3. 

The further structure of our paper is as follows. In Section 2, we present and discuss the 
main results. In Section 3, the averaging transformation is used to simplify system (1.1). The 
necessary properties of the unperturbed system (1.4) are discussed in Section 4. In Section 
5, we introduce local coordinates in the vicinity of the unperturbed invariant torus. Further, 
we show the existence of the perturbed integral manifold and study the dynamics on this 
manifold in Section 6. Remaining proofs of main results are given in Section 7. Finally, in 
Sections 8 and 9 we illustrate our theory by two examples. 

2. Main results. In the co-rotating coordinates x{t) = e^"'^^y{(5ot) the unperturbed equa- 
tion (1.4) reads as 

/3o^ = f{y) - aoAy. (2.1) 

The quasiperiodic solution (1.5) to (1.4) now appears as 27r-periodic solution y{^p) = xo{ip) to 
(2.1). The corresponding variational equation around this solution is 

f!^o^ = {f'ixom-aoA)y. (2.2) 

It is easy to verify (see Section 4), that (2.2) has two linear independent (because of assumption 
(1.6)) periodic solutions 

9i(V') := ^(V'), q2W ■■= Axoii^), (2.3) 

which correspond to the two trivial Floquet multipliers 1 of the periodic solution xq to (2.1). 
One of these Floquet multipliers appears because of the 5^-equivariance of (2.1), and the other 
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one because (2.1) is autonomous. Prom the exponential orbital stability of (1.5) it follows that 
the trivial Floquet multiplier 1 of the periodic solution xq to (2.1) has multiplicity two, and the 
absolute values of all other multipliers are less than one. Therefore, there exist two solutions 
Pi{ip) and P2(V') to the adjoint variational equation 

Poj^ = -{nxo{i;)f + aoA)z (2.4) 

with 

pJWqk{^) = Sjk (2.5) 

for all j,k = 1,2 and "0 G l^j where 6jk = 1 for j = k and 5jk = otherwise. 

2.1. Non-degenerate case. Using notation (1.9) we define a 27r-periodic function Gi : 

M M by 

27r 

GiW -.= ^1 pfii; + e)giixoi^ + e),9)d9 (2.6) 


and its maximum and minimum as 

Gt ■= max Gi(ih), G7 := min Gi(ih). 

^ t/.e[0,27r] ^ VG[0,27r] 

For the sake of simplicity we will suppose that all singular points of Gi are non-degenerate, 
i.e. 

G'Kip) / for ah with G;(0) = 0. 

This implies that the set of singular points of Gi consists of an even number 2N of different 
points. The set of singular values of Gi will be denoted by 

5i := {Gi(V'i), . . . , Gi{^2n)}, where {V^i, . . . , V27v} = {V' e [0, 27r) : G[{^P) = 0}. 

Our first result describes the behavior (under the perturbation of the forcing term in 
(1.1)) of Ti X M, which is an integral manifold to (1.4), as well as the dynamics of (1.1) on the 
perturbed integral manifold to the leading order. 

Theorem 2.1. For all /3i < P2 there exist positive 7* and a^, such that for all a,/3, and 7 
satisfying 

a>a,, /3i</3</32, < 7 < 7,, (2.7) 

system (1-1) has a three-dimensional integral manifold 9Jli(a, /3, 7), which can he parameterized 
by tp,ip,t £ M. in the form 



X = e^'^xo(V') + 7e^'^f^i (V-, 7) + -U2 (V', l^t, at, 7, 0"^ , (2.8) 

a 

where functions Ui and U2 are C'^^ smooth and 2-7T-periodic in if), ip, (it, and at. 
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The manifold 97li(a,/3,7) is orhitally asymptotically stable and solutions on this manifold 
are determined by the following system 

^ = /3o + 7pJ{iP)gi{xo{ij),m + 7^Su + -S12, (2.9) 
at a 

^ = ao + lplWgi{xoW,/3t) + j^S2i + -S22, (2.10) 
at a 

where functions Sjk = Sjk{ip,^, f3t,at,'y,a~^) are C'^^ smooth and 27r-periodic in il),ip,j3t, 
and at. 

In addition, for any e > there exist positive ai,7i, and 6 such that for all a, j3, and 7 
satisfying 

a > ai, < 7 < 71, 7Gr < /3 - /3o < 7^^ dist ^i) > e (2.11) 

the following statements hold: 

1. System (1.1) has integral manifold 9Jti(a,/3,7) of the form (2.8) and even number of 
two-dimensional integral submanifolds 9Tj C 5[Ri(a, /?, 7), j = 1,...,2A^, N = N(a,P,^) < N, 
which are parameterized by {p,t £ M, in the form 

X = e^'^xoi^t + 'dj) +jVij{^,f3t,at,-f,a-^) + ^V2j{(p, (3t,at,j,a~^). (2.12) 
The dynamic on the manifold 9Tj is determined by a system of the type 

^ = ao + jWj {(f, pt, at, 7, a^^) . 

Here dj are constants and functions Vij, V2j, and Wj are C'^^ smooth and 2t: -periodic in 
ip, (p, (3t, and at. 

2. Every solution x{t) of system (1.1) that at a certain moment of time to belongs to a 
5 -neighborhood of the torus T2 tends to one of the manifolds 9Tj as t ^ +00. 

One of the main statements of this theorem is that under conditions (2.11), within the 
perturbed manifold 9Jti(a, /3, 7), there exist lower-dimensional "resonant" manifolds cor- 
responding to the frequency locking. These manifolds attract all solutions from the neighbor- 
hood of T2. Hence, the asymptotic behavior of solutions is described by (2.12) and has the 
modulation frequency (3. One can observe also from (2.12), that the perturbed dynamics is, 
in the leading order, a motion along Ti with the new modulation frequency. The following 
theorem describes these locking properties more precisely. 

Theorem 2.2. For any e > and ei > there exist positive 71,01, and 5 such that for all 
parameters (a,/?, 7) satisfying the conditions (2.11) and for any solution x{t) of system (1.1) 
such that dist(x(to), Ti) < 5 for certain t^ & there exist cr, T G M such that 

inf \\x{t) - e^'^xo(/3t + < ei for all t > T. (2.13) 
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Generic case 




Figure 2.1. Cross-section of the locking regions with respect to the frequency /3 and amplitude 7 by a = const 
of the external perturbation for the generic case gi ^ 0. 



The conditions for the locking (2.11) depend on the properties of the function Gi. For 
the case that Si consists of two numbers Gj~ and Gf only, i.e. that Gi has only two singular 
values, the set of parameters satisfying these conditions is illustrated in Fig. 2.1 (left) for a 
fixed value of a > ai . The admissible values of the parameters (3 and 7 belong to a cone with 
linear boundaries 

/3 = /3o + 7(G'r+e) and/3 = /3o + 7(G^-e), 

bounded from above by 7 < 71. 

If Gi has more than two singular values, then the corresponding regions in the parameter 
space 

|/3-/3o-7Gi(Vi)| <7e 

should be excluded. Here ^pj are the additional singular points of Gi. An example is shown 
in Fig. 2.1 (right) for the case with two additional singular points. 

2.2. Degenerate case. In this subsection we suppose that 

(2.14) 



1 

— / g{x, tp, ip)dip = for ah x G M*^ and G 
2vr Jo 



In this case the functions gi and Gi, which are defined in (1.9) and (2.6), are identically zero 
and, hence, cannot give any information about locking behavior like in Theorem 2.1. Instead, 
the following functions g2 xM. ^ R" and ^2 : M — >• M will define the locking: 



and 



27r 



27r 



2-iT 



dui 
dx 



{x,'ip,(p)g{x,i;,ip)d(p 



G2W / pUV' + %2(2;o(V' + 



(2.15) 
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Here 

g{x,ip,e)de 

is the antiderivative of g{x, ip, ip) satisfying 

r-2n 

/ ui{x, tj), ip)dip = 0. 
Jo 

Now we proceed as in the non-degenerate case: We denote 



Gt := max G2(it'), Go := min G2(ip) 

^ V6[0,27r] ^ V6[0,27r] 

and suppose that all singular points of G2 are non-degenerate 

G'2{i)) / for all ^ with G'^i:^) = 0. 

This implies that the set of singular points of G2 consists of an even number 2N of different 
points. The set of singular values of G2 will be denoted by 

52 := {G2{iJi), . . . , G2(V'2iv)}, where {V'l, • • • , 4'2n} = {V' e [0, 27r) : G'^{iP) = 0}. 

Our next result describes the locking behavior of the dynamics close to Ti in the degenerate 
case. 

Theorem 2.3. Suppose (2.14) holds. Then for all /3i < (32 there exist positive constants /i* 
and a^, such that for all a, (3, and 7 satisfying 

/3i < /3 < /32, 7Va < /^*, « > a*, (2-16) 

system (1.1) has a three-dimensional integral manifold 9Jt2(«) /3, 7), which can he parameterized 
by ip,ip,t G M in the form 

X = e^'^xoW + ^?7i fv,</',/3t,at,— (i^, if, ^t, at,^, -) , (2.17) 



a \ a a J a \ a a 



where functions Ui and U2 CLre G^^^ smooth and 2TT-periodic in ip, ip, fit, and at. 

The manifold dJl2 {a, (3, is orbitally asymptotically stable and dynamics on this manifold 
is determined by the following system 

^ = /3o + ^pJ{iP)g2{xoW,m + ^^11 + -^812 + :^5i3, (2.18) 



2_ 


S\\ 


+ 




5'l2 


+ 


7^ 


S^ 








oP- 










7 


S21 


+ 


1' 


'S'22 


+ 


1' 




a2 






a2 






a2 





§ = ao + ^pl{^)g2{xoi^),l3t) + ^^21 + ^^22 + ^^23, (2.19) 
at a 

where functions S^k = ("01 V') /^i, 7^/ct, 1/a) are C^~^ smooth and 2n-periodic in ■ip,(p,/3t, and 
at. 
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In addition, for any e > 0, there exists some positive ai,ci, C2, and 5 such that for all 
a, 13, and 7, satisfying 

a>ai, - < 7 < C2^/^, —G^ < /3 - /3o < —G+, dist ( ^{f3 - f3o), S2) > e, (2.20) 
a a a \7^ / 

then the following statements hold: 

1. System (1.1) has integral manifold 9Jt2(a,/3,7) of the form (2.17) and even number of 
two-dimensional integral submanifolds £j C j = 1,...,2N, N = N{a,P,^) < N, which 
are parametrized by (p,t G M in the form 

X = e^'^xoiPt + i^j) + -^Vij + ^V2j + -Vsj + — F4i, (2.21) 
a a a 

where C'~^ smooth functions Vij = Vij{(p, j5t,at,^'^ / a^\/ a^\/ a^y) are 2'K-periodic in (p,(3t, 
and at. The system on the manifold reduces to 

^ = ao + ^Wji + ^Wj2 + ^Wj3, 
at a a-'' 

where functions Wji = Wji{ip, I3t,at,j'^ /a,l/a,l/a^) are C'~^ smooth and 2iT-periodic in 
if, f3t, and at. 

2. Every solution x{t) of system (1.1) that at a certain moment of time to belongs to 
5 -neighborhood of the torus T2 tends to one of the manifolds 2,j as t ^ +00. 

Theorem 2.3 gives verifiable conditions on the parameters (a, /3,7), for which the locking 
of modulation frequency to the modulation frequency /3 of the perturbation takes place for 
the case when the degeneracy condition (71 = is fulfilled. These conditions are given by 
(2.20) and differ from the conditions of locking in the generic case given by Theorem 2.1. 
The locking phenomenon in the leading order looks similarly in both cases, i.e. the solutions 
tend asymptotically to Ti and modulation frequency (3. More precisely, the following theorem 
holds. 

Theorem 2.4. For any e > and ei > there exist positive ai, ci, C2, and 6 such that for 
all parameters (a,/3,7) satisfying the conditions (2.20) and for any solution x{t) of system 
(1.1) such that dist(x(to), Ti) < <J for certain t^ gM. there exist o", T G M such that 

inf ||x(t) -e^^xo(/3t + a)|| < ei for all t > T. (2.22) 



Let us illustrate the set of parameters (2.20) leading to the locking and compare it with 
the generic case. For a fixed large enough a, the region in the parameter space (/3, 7) has 
again a shape of a cone as in Fig. 2.1, but now the cone has tangent boundaries 




(2.23) 



leading to a smaller synchronization domain for moderate values of 7, see Fig. 2.2. Small 
values of 7 are even excluded 7 > ci/a. But, on the other hand, the synchronization is now 
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Degenerate case gi = 







a 



a 



Figure 2.2. Cross-sections of the locking regions with respect to the frequency (3 and amplitude 7 by 
a — const of the external perturbation for the degenerate case gi = 0. 

allowed for large values of 7 up to C2\/a. This means practically, that the locking occurs not 
only for small but also for large amplitude perturbations, contrary to the case gi ^ 0. 

Again, when the singular set 5*2 of G2 contains more than 2 points, then the corresponding 
regions in the parameter space 



should be excluded. Here V'j are additional points from 82- Example in Fig. 2.2 (right), shows 
the case with two additional singular points. 

Finally note that the locking phenomenon described in [19] for a simple model with sym- 
metry corresponds to the case when an additional degeneracy takes place with 52 = 0. As it 
is shown in [19], the locking cone becomes even more high (7 < /U*a) and thin with the slope 
proportional to 7^/a^ instead of 7^/a in (2.23). 

3. Averaging. We perform changes of variables which depend on the average of the per- 
turbation function g with respect to fast oscillation argument at. As the result of these trans- 
formations we obtain an equivalent system, where the fast oscillation terms have the order of 
magnitude of 7/a and 7^ /a and smaller. The principles and details of the averaging procedure 
can be found e.g. in [2, 10, 21]. This transformation has the form 



/3 - /3o - ^i^j 




1' 
— e 
a 



X = 



xi -\ — ^1(2;, /3t, at). 




a 



Inserting (3.1) into (1.1), we obtain 



dxi 7 / dui dx dui 




+ 7 



dui 



f [xi + -ui) + jg{x, l3t, at) 



dt a\dx dt d{(3t) 



d{at) 



a / 
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Accordingly to the idea of the averaging method, the terms of order 7 depending on the high 
frequency at should vanish due to the choice of ui . This leads to the condition 

g{x,l3t,at)-gi{x,/3t), (3.2) 



diat) 
where 



1 

gi{x,(Jt) = ■ 

Hence 



2-K 



5(i(x,/3t) = — g{x,l3t,(p)dip. 



/•at 

ui{x,/3t,at) = / {g{x,/3t,ip) - gi{x,/3t))dip 

J tfo 

is a periodic function in at and /3t. Strictly speaking, the above integral is considered com- 
ponentwise, and the initial points (po can be different for each component of vector-function 
Ui. Additionally, we choose the initial points (poj in such a way that 

ui{x, (3t, ip)dip = 0. 



The resulting averaged system reads 

dxi 'y 0'^ 

— = /(xi) + 75i(xi,/3t) + ^Ai{xi,f3t,at) + ^A2{xi,(3t,at) 
dt a a 

2 3 

+2-Xi{xi,l3t,at,^) + 2-X2{xi, /3t,at,^), (3.3) 
a a^ a 



where 

MxuPt,at) = -u,-—f-g^^f^, 

A / n \ 9ui dgi 
A2{xi,pt,at) = --Q^9+ "^^1- 

Functions Aj{xi,j3t, at) and Xj{xi,f3t, at, ^), j = 1, 2, are C'~"'^-smooth and 27r-periodic in (3t 
and at. Here we use the fact that for small 7/a the change of variables (3.1) is equivalent to 

X = xi + —ui{xi, I3t,at,—) (3.4) 
a a 

with smooth bounded and periodic with respect to /3t and at function ui. 

Note that the averaged term gi{x,'ip) is equivariant with respect to the group action e^^ : 

gi (e^«x, ^) = i- ^ g{e^^x, tp, ip)d^ (3.5) 
= ^ r e^^g{x,il^,yp - Od^ = e^^gi{x,i;). (3.6) 
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If gi{x, -i/;) = we perform the second averaging change of variables 

7 7^ 
Xi = X2-\ 7tU2i{xi, /3t,at) H :tU22{xi, f3t, at). (3.7) 

The functions U2i and U22 are selected from the conditions that the terms of the orders 7/a 
and 7^/a, which depend on high frequency at, vanish. This leads to 

—-— = Ai{xul3t,at), —— = A2(xi,/3t,at) - g2{xi,l3t), (3.8) 
o[at) o[at) 

and 

52(xi,/3t) = ^ j^J A2{x^,pt,^)d^ = -i-£"M^^^5(xi,/3t,<^)d(^. (3.9) 

We have used in (3.8) the following property of Ai: 

1 /■^'^ 

— Ai{xi,pt,ip)dip = 0. 

The second averaging function g2ix,Pt) is also equivariant with respect to the symmetry 
group action 

g2{e^^x,'ilj) = e^^g2{x,Tp). 
This can be seen from the following calculations 

<72(e^«x,V) = -^1^^ (£^1!^''^^'^'^'^^'^^) 9{e^^^^^^^)'^V 

c2tt / Pip 



-e 



After the second averaging, in the case gi{x,f3t) = , the system admits the form 



dx2 "7^ 7 7^ 1 7*^ 7^ 1 

-TT = f{x2) H g2{x2,Pt) + ^g3(2;2, — , -) + -^g4{x2, ^t, at, —, - 

dt a a a a'^ a a 
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y (4.2) 



where C'~^-smooth functions 33 and are 27r-periodic in at and j3t and function g2{x2, f^t) 
is 27r-periodic in pt and equivariant with respect to action e^^. Note that by (3.1) and (3.7) 
the variable x is expressed by X2 in a form 

X = X2 + —U2{x2, (3t,at,—,—) (3-11) 
a a a 

with smooth bounded and periodic with respect to f3t and at function U2. 

4. Useful properties of the unperturbed system. In this section, we consider useful 
properties of the linearized unperturbed system and introduce an appropriate basis (matrix 
$0)) which locally splits the coordinates along the invariant torus (1.7) and transverse to it. 
Further, this basis will be used in section 5 for the introduction of appropriate local coordinate 
system. 

Since the unperturbed system (1.4) has quasiperiodic solution x{t) = e^'^°^X(){P()t) then 

^^^M + aoAxoiM = fixoiM)- (4.1) 
The corresponding variational system 

dy ^ df{x) 
dt dx 

has two quasiperiodic solutions 

Ae^-»*xo(/3ot) and e^""*^^^^. (4.3) 
The following properties of the Jacobian follow from the equivariance conditions 

dm ^ ^-Ae dfie^'-) ^Ae^ (4.4) 
dx dx ' 

Af{x) = ^^Ax. (4.5) 
The latter conditions and the change of variables y{t) = e^°"^^w{l3Qt) in (4.2) lead to 

where ■0 = Pot- The linear periodic system (4.6) has two periodic solutions 

Qi(V') = ^^, q2W=AxoW, (4.7) 

see (4.3). 

Let 3 be the trivial vector bundle 3 = (IR" x 'ri,Ti,p), where p is the natural projection 
onto Ti. Consider corresponding to (4.6) linear skew-product flow with time ip and xq E M", 

7r(V', xo,'il^o) = (J7(V', i^o)xo, + -00) , (4.8) 
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where r2('0,'i/'o) is the fundamental solution of (4.6) such that r2('0o,V'o) = I- The vector 
bundle 3 is the sum of two sub-bundles 3i and 32, which are invariant with respect to the 
flow (4.8). The two-dimensional bundle 3i consists of periodic solutions of (4.6) spanned 
by two linearly independent periodic solutions (4.7). The solutions from the complementary 
bundle 32 tend exponentially to zero as t — t- oo. Since the bundle 3i is trivial, the bundle 32 is 
stably trivial. By [11, p. 117], any stably trivial vector bundle whose fiber dimension exceeds 
its base dimension is trivial. Therefore, if n > 3, the {n — 2)-dimensional bundle 32 is trivial 
and there exists a smooth map <I>o ^ I"! ~^ /^(M""^, M"), which is isomorphism between 32 and 
Ti X M"-^ In the case n = 3, is a 3— dimensional vector function, whose existence can 

be shown by direct analysis. 

By construction, n x n-matrix 

^i(V') = {gi(^),'?2(^),$o(^)} (4.9) 

is non-degenerate for all ij:. By the change of variables w = ^i{tjj)z system (4.6) is transformed 
to system 

^ = diag|o,0,^L, (4.10) 



where (n — 2) x (n — 2) matrix H{-ip) is 27r-periodic in and subsystem 

is exponentially stable. 

Since linear periodic system (4.6) has two nonzero linearly independent periodic solutions, 
the adjoint system 

^ = -BiWw (4.11) 

has also two nonzero linearly independent periodic solutions Pio(V') and P2o(^) [6]- Corre- 
sponding to (4.11) linear skew-product flow on 3 = (K*^ x Ti,Ti,p) has two invariant sub- 
bundles 3i and 32- Two-dimensional bundle 3i consists of periodic solutions of (4.11) and 
solutions from 32 tend exponentially to zero as t — )■ — oo. As in the case of system (4.6), 
for the adjoint system (4.10) there exists a smooth map $oo '■ Ti — £(]R"~^, R"), which is 
isomorphism between 32 and Ti x IR""^. Therefore, the matrix 

^i(^) = {pwW,P2oW,^ooW}, 

is non-degenerate for all ip. 

Let us show that the following holds 

l>f(V')a>i(^) = diag{Ci,C72(^)}, (4.12) 

where Ci is a constant non-degenerate 2x2 matrix and C2{ip) is a (n — 2) x (n — 2) non- 
degenerate periodic matrix. Indeed, the scalar product of any solution of (4.6) and any 
solution of its adjoint system is constant for all t/j [6]. Hence, pfQ{'4^)<lj{ip) = const. Moreover, 
the scalar product {ij))w{^lj) of any solution w{ijj) from 3i of system (4.6) and a solution 
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■w{ip) from 32 of its adjoint system is zero for all ip, since w{'ip) is periodic and ■w{'ip) — )• for 
^lJ — )• — oo. Therefore, for all ^, the fibers 3i(V') ^ind 32 ("0) over point ip are orthogonal as 
subspaces of M". Similarly, the subspaces 3i(V') ^^'i 32 (V') ^-re orthogonal. This leads to the 
block-diagonal form of the product (4.12). 

From the equality (4.12) follows that ^fHV') = diag {Cf \ C2"^(V')} ^H?/')- Therefore, 
the matrix <I>j~^(^) has the form 

'^I'W = {pi(V),P2(V'),^'o(V')}^, (4.13) 

where pi{ip) and P2(^) are periodic solutions of (4.11) such that pj {'ip)qi{'4^) = p'2 {'^)Q2{ip) = 1 
and pj {ip)q2{'4^) = Pi"(^)92(V') = for all ip. 

Taking into account (4.10), it can be verified that 

^'^^(^)+$,(^)diag|o,0,^Ui^o(^)'^i(^). 



dip { /3o 

Then the n x {n — 2)-matrix <^>o(^) satisfies the relation 

+ l$o(^)i^(V') = Boimoii^). (4.14) 

dip po 

5. Local coordinates. In this section we write systems (3.3) and (3.10), which appear 
after the averaging transformation, in the local coordinates in the vicinity of the invariant 
torus. Systems (3.3) and (3.10) have the form 

— = /(j;) + fi'^g{x,/3t) + fi'^e'^ri{x, I3t,at, fi,e) + fj.e^r2{x, /3t,at, fi,£), (5.1) 

where / and g are 5^— equivariant, i.e. f{e^'^x) = e^^ f{x) and g{e^'^x, fit) = e^^g{x^ fit). In 
particular, system (3.3) can be obtained from (5.1) by setting 7 = /i^, 1/a = e^, ^ = (71, and 
r2 = 0. System (3.10) has the form (5.1) if = /i^, 1/a = e^, and g = 52, respectively. Here 
the perturbation terms fi'^e^ri and fie^r2 are defined from (3.3) and (3.10) in a straightforward 
way. In this and the following section, we consider n and e as independent parameters. 

We introduce new coordinates ip, f and h instead of x in the neighborhood of the two- 
dimensional torus Ti by the formula 



X = e 



''^{xoi^) + Mi^)h), (5.2) 



where h S ^, \\h\\ < ho and the matrix ^o{iJj) was defined by the trivialization of bundle 
32, see (4.9). 

Substituting (5.2) into (5.1), we obtain 

= /(e^^(xo(V') + Mi^m + ^i2^(e^^(a;o(V') + <J>o(V')/i), /?*) 
+ /i2eVi(e^'^(xo(V') + $o(V')^), /3i, ^, e) 

+ /ieV2(e^^(xo(V') + Mi^)h),f3t, at, /x, e). (5.3) 
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The following transformations allow to split the variables h, ip, and (p. Let us write (5.3) in 
the form 



\ di; dil: )\dt ^° 



dip 



+ A (xo(V) + M^)h) 37 - ao + «'o(V') -n - H{i,)h 



dt 



dh 



dt 



+ /i2e2e-^'^ri(e^'^(xo(^) + Mm,l3t, at, fi, e) 

+ fie^e-'^'^r2ie^'^ixo{i^) + $o('0)/i), at, ^, e). (5.4) 
Taking into account (4.1) and (4.14) we have 



f -/3o 

f - ^(V^)/i 



/(xo + 'foh) - f{xo) - ^^^o/i 



+ ^'^g{xo,(3t) + fi'^g{h, V', /3t)/i + iJ^e'^rs{h, V', 95, /3t, at, /i, e) + ^e'^r4^{h, ip, 93, /3t, at, ^, e). 



(5.5) 



where 



^i(^) 



dtp 

rf^o(^) 



AxoiiP) $o(V') 



g{h, V, /3t)/i = 5(xo(V') + Mi^)h, f3t) - g{xoi^), Pt). 

Since by our construction det$i('i/') 7^ for all tj: £ Ti, the matrix ^i{ip) + <I>2(V'',^) is 
invertible for sufficiently small h and 



Hi{h,tP) 



($i(V') + $2(V',/i))- =^1 (^) + ^(^>V')= PiW + H2{h,tp) , 



(V^) 



i/3(/l,V') 



where C' ^-smooth function H[h,ip) = 0{\\h\\) is periodic in tp. 

Therefore, system (5.5) can be solved with respect to derivatives dip/dt, dip/dt and dh/dt : 

^ = /5o + fi^pf{tp)g{xo,Pt) + Ru + fJ^^Ri2 + fJ^^e'^Ru + f^e^Ru, (5.6) 
^ = ao + fj.^p^{tp)g{xo,l3t) + R21 + ii^R22 + Ai^e^^23 + ^e^i?24, (5.7) 



dt 



(5.8) 
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where 

Rji = Rjiih,ij,fi) = {pJii^) + Hjih,i;))Fiih,^), 

Rj2 = Rj2{h, V, f-i) = (pJ(V') + Hj{h, xP))g{h, V', Pt)h + Hj{h, iP)g{xo, /St), 

Rjk = Rjk{h,ip,ip,f3t,at,fi,e) = (j>J{ip) + Hj{h,ip))rk{h,'ip,ip, f3t,at, fi,e), 

Rm = Rm{h,^,li) = {^oW + H3{h,i;))Fi{h,i;), 

Rm = Rm{K ^, = ('^o(V') + H^{K ^))g{xo{'^) + M^)h, Pt), 

Rok = Rosih, ■0, ^, /3t, at, fj,, e) = ($o(V') + ^3(^, '4'))rk{h, tp, cp, (3t, at, e), 

F,{h, ^) = f{xo + - f{xo) - ^^^oh = 0{\\hf) 

with j = 1,2, k = 3,4. 

The foUowmg lemma estabhshes the existence of the pertm'bed manifold. Note that the 
similar lemma has been proved in [19] for the system with constant matrix H and somewhat 
different parameter dependences. 

Lemma 5.1. For fi G (0,/io) one? e £ (0,eo) with sufficiently small fiQ and Eq, the system 
(5.6) - (5.8) has an integral manifold 

mf,,e = {{h,i^,ip,t) : h = u{i;,ip,^t,at,fi,e), (V,<^) GTi xTi,t gM}, (5.9) 

where the function u has the form 

u{ip, (f, I3t, at, fi, e) = fi'^uo{il), /St, fi) 

+ e^fiui{il;,ip, (it,at, fi,e) + e'^ fi'^U2{ip,(p, (3t,at, fj,,e). (5.10) 

Here functions uq,ui and U2 are C^~^-smooth, lir-periodic in 'ip,ip,f3t and at, satisfying 
\\uj\\(ji-3 < Ml, j = 0,1,2, with positive constant Mi, which does not depend on /i, e, and 
a > aQ. Here \\.\\ci-3 is the norm of functions from C'^^(T4) with fixed parameters /i and e. 

The integral manifold 9?t^^e is asymptotically stable, i.e. there exists vq = VQ{jJL,£Q) such 
that for every initial value {h,ip,(p) at time r with \\h\\ < vq there exists a unique {ij)Q,(pQ) 
such that 

\\N{t, T, h, ip, (f) - N{t, T, n(V'o, V^o, /3t, ar, /i, e) , ipo , ipo)\\ 

< Le-^^^-^'^\\{h,^,ip) - (u(0o,<^o,/3r,aT,/i,e),V'o,V5o)||, (5.11) 

where t > t, L and k are some positive constants not depending on a, e, N{t, t, h, ip, ip) is 
the solution of the system (5.6) - (5.8) with an initial value N{T,T,h,%lj,ip) = (h,'ip,ip). 

Proof. By introducing new variables (" = (p, (i, (^2)7 Ci = f^t^ C2 = at, and new parame- 
ters A = {rii,r]2,r]^, fJ., e), r]i = y?, r/2 = e^Ai^, ^3 = e^/i in system (5.6) - (5.8), we obtain the 
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following autonomous system 



^=i/(V')/i + Qo(/i,C,A), (5.12) 
^ = L^o + Q(/i,C,A), (5.13) 



where = {(3o,ao, f3,a) and Q = (Qi, Q2, 0, 0), 

Qo = Roi + viRm + mRo3 + %-Ro4 

The corresponding reduced system has the form 



Qj = mpjg + Rji + 'niRj2 + mRj3 + ??3^i4, j = i, 2. 



§ = «(,>)^ |=.„. (5.14) 



By construction, this system is exponentially stable, i.e. for all ip ^ Ti the fundamental 

dh 
dt 



solution Clo{t,T,ip) of system ^ = H{j3Qt + il})h satisfies 



||^^o(*,-r,V')|| < re-'^o^*-^), t > T, (5.15) 

where £ > 1 and kq > 0. Note that matrix H{ijj) and the fundamental solution il.Q{t,T,'ip) 
depend only on the first coordinate of vector C = {ip, 93, Ci, ^2)- 

Denote Ixq = {A : ||A|| < Aq}. Let J^p be the space of Lipschitz continuous functions 
w : T4 X I\q — )• R" such that \\w\\c < p, Lip^ w < p for all A G /aq, where Lip^tf is the 
Lipschitz constant of w with respect to the first argument C- 

In order to prove the existence of the invariant manifold for (5.12) - (5.13), we consider 
the mapping T : Fp ^ Fp 

T{W){C,X)= f n{T,OQo{w{Cr,X)Xr,X)dT, 
J —00 

where Cr = {iPt,^tXit,C2t) is the solution of (5.13) for h = 'w{C,X) with initial condition (, 
i.e. 

^=u;o + Q{w{C,X),C,X). 
r2(f, C) is the fundamental solution of system 

§ = H{^ljt)h, (5.16) 

with 17(0, C) = I- As follows from [28] (Lemma 6.3), system (5.16) is exponentially stable as 
a small perturbation of system (5.14), therefore there exist Aq > and po > such that for 
A G Ixo and w G Tp^ 

||f](t,C)|| </:ie-"^*, t>0, (5.17) 

with Ci > 1 and < ki < kq. 
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Analogously to [19] we prove that the map T(w) : Fp ^ Fp \s & contraction for small 
enough pQ and Aq- Hence, the mapping T{w) has unique fixed point wo{C,,X). 

For proving C'~^ smoothness of integral manifold wq{C,,X) we use the fiber contraction 
theorem [25], [3, p. 127]. In particular, for the proof of C^-smoothness with respect to C,, we 
introduce the set J-"^ of all bounded continuous functions $ that map T4 x Ix(^ into the set of 
all n X 4 matrices. Let J-p denote the closed ball in with radius p. For w ^ Fp, we consider 
the map T^{w, <1>) : Fp x Fj, ^ Fp, which is defined as follows 

+ ^q°<'"K^;^^)- <-'> MC^,X)W^ + ^>^Kr, A)) dr. (5.18) 

where Ct and Wt are solutions of the system 

^=^o + QMC,A),C,A), (5.19) 

dw ^ dQ{w{cxU,x) ^ ^ aQ(^CA),c,A) ^ _ 

dt dC dh ' 

Analogously to [3, 19] we prove that the map 

{w,<^) ^ {T{w),T^{w,^)) (5.21) 

is continuous with respect to w and $ and is a fiber contraction. It has unique fixed point 
{wq,wi) which is globally attracting and bounded uniformly with respect to a £ [ao,oo). 
Repeating [3, p. 296] one can show that wq is continuously differentiable and Dwq = wi. 

The smoothness up to C'~^ can be improved inductively. The continuous differentiability 
with respect to A is proved analogously. 

Since the invariant manifold h = w^C, A) for f?i = //2 = % = equals to zero /i = 0, it can 
be represented as 

h = r]iWo{i/j,Ci,f^) + mwi{^p,ip,Ci,C2,>^) + r/3W2(^/^,(/?,Ci,C2, A). 

Note that wq does not depend on (2, f/2i%5 and e, since system (5.12) - (5.13) is independent 
on (^2 for ?72 = r/3 = e = 0. Taking into account the definitions of r/i, r]2, and Cii C21 we 
obtain that the invariant manifold of (5.6) - (5.8) has the form (5.10). 

The remaining proof of (5.11) is analogous to the proof of Lemma 5.1 in [19]. ■ 
Remark l.If r2 = in system (5.1) then i?04 = 0, Ru = R24 = in system (5.6) - 
(5.8) and ui = in the expression for the integral manifold (5.10). This corresponds to the 
non-degenerate case described by (3.3) with non- vanishing averaged term gi{xi,f3t). 
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6. Investigation of the system on the manifold. Substituting the expression for the 
invariant manifold (5.10) into equations (5.6) - (5.8), we obtain the system on the manifold 

^ = /3o + /i^pf (^)^(xo(V'), /3t) + /Sii(V', f3t, fi) + 

+£^H^Si2{'>p, ^, (3t, at, fi, e) + e^fiSi3{ip, ip, (3t, at, /i, e), (6.1) 

IT 4 
ttO + Ai P2 {'^)g{xQ{llj),l3t) + /X 52l(?/', Pt, fl) + 



dt 

+e^/i^5'22(V'i ^, Pt, at, /i, e) + e^iJ,S23{ip, ^, Pt, at, /i, e), (6.2) 

where C'~^-smooth functions Sji are 27r-periodic in tj), ip, fit, and at. 
Now we will use the closeness of the frequencies /3o and /3 

/3 - /3o = /i'A 

and obtain conditions for the locking of the variable if) to the external frequency (3. For this, 
we introduce new variable "01 in (6.1) - (6.2) according to the formula 

ij) = /3t + ipi. 

In the obtained system 

^ = -fi^A + /iV(/3t + ^i)g(xo(/3t + i>i),f3t) + /5ii(/3t + iPi,pt, ii) + 

+e^f/Si2{/3t + V'l, 9?, l3t, at, fi, e) + e^fiSisO^t + ^fl,lp, j3t, at, fx, e), (6.3) 
^ = ao + ^^P^(/3t + ifiMMl^t + Vi), /3t) + /S2i(/3t + ^1, /3t, ^) + 

+e^/i^S'22(/3t + V'l, 9?, ai, fJ-, e) + e^/^5'23(/5i + 'i/'i) V', a*, A*) e) (6-4) 

the variable ifi is slow and one can perform the following averaging transformation with 
respect to /3t 

,,2 



^1 = ^2 + J I [pi{i + ^i)-g{xo{i + ^i),0-G{iJi)]di, 

^2 rPt 



where 



^ = ^2 + J I [P2(e + ^i)9(^o(e + v^i),6-^2i(V'i)]d?, 



G(Vi) = ^ I^^P^^^ + V'i)5(xo(^ + 



1 



27r 



52i(V'i) = ^ / P2{C + i^i)9{xo{^ + i^i),OdC- 
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After this transformation system (6.3) - (4.13) takes the form 

^ = -A^2 ^ ^2^(^^) ^ /5n(V2, ^i) (6.5) 

+e'^fi'^Si2{ip2, ^2,/3t, at, fi, e) + e^iJ,Si3{tp2, ^2,Pt, at, /i, e), 

^ = ao + iJ?S2i{i)2) + /i^52i(?/'2, Ai) (6.6) 
at 

+e^/i^5'22('i/'2, ¥^2, ai, A^, e) + e^/^5'23(V'2, ¥52, ai, A*, e), (6.7) 

where the functions in the right hand side are C'~'^-smooth and periodic in iIj2, (p2, f3t, and at. 
The corresponding averaged system allows splitting off of the dynamics for 1^2 variable 

^ = /i2(-A + G(^2)), (6.8) 

which can be described completely. In particular, it has equilibriums, which are solutions of 
the equation 

A = G(0. (6.9) 

Such equilibriums exist if 

A = ^G[G-,G+], 

where G~ and are minimum and maximum of the periodic function G{(,), respectively. 

If A € [G~ , G"*"] is a regular value of the map G, then there exist even number of solutions 
? = j = 1, . . . ,27V, of equation (6.9) such that G'('(?°) / 0. In general, N depends on A. 
Since the signs of every two sequential values G'(?9°) and G'{■l^^_^_l) are opposite, the half of 
these equilibriums are stable and another half is unstable. Every such equilibrium corresponds 
to an integral manifold of (6.5)-(6.6) with the same stability properties. 

Lemma 6.1. Assume A G [G~,G'^] such that (6.9) has 2N nondegenerate solutions ^ = t?^, 
j = 1, . . . , 2N . Then there exist /io > and cq > such that for allO < fi < fiQ and e < cofi^^^ 
system (6.5) - (6.6) has 2N integral manifolds 

= {(^2,(^2, i) : i^2=^^j+Vj{m,^t,at,fi,e), ^2&Ti,t€R} , (6.10) 

where 



^3 

Vj = fi'^Vjo{f3t, fx) + e^Vji{ip2,Pt,at,^i,e) H Vj2{f2,^t,at,fi,e), 



with G'"^ smooth, periodic in ip2,(3t,at functions vjk, such that \\vjk\\(ji-3 < with the 
constant M3 independent on fj,,e, and a > a^. 

The manifolds Il2k, k = 1,...,N, are exponentially stable in the following sense: there 
exists 60 such that if \xp20 — i?2fcl — '^0 o^c? tpo ^ Ti, then there exists a unique (poi such that 
for t > to the following inequality holds 

\'>p2{t,to,Tp20,^o) -i^2{t,to, ??2fc + V2k{Voi,Pto, ato, fl, e),ipoi)\ 

+ \V2{t,to,Tp20,^o) - ip2{t,to,i}2k +V2k{Voi,Pto,ato,fl,e),(poi)\ 

< /:2e-'^''=2(*-*°) {\ipo - (^oil + IV'20 - ^2k - V2kiiPoi,/3to,ato,ti,e)\) , (6.11) 
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where C2 > 1 and K2 > is independent on a, /U, and e. 

The manifolds n2j!c_i, /c = 1, N ^ are exponentially unstable in the following sense: there 
exists 5q such that if \tp20 — '^'2k-i\ — ^0 '^'^'^ G Ti, then there exists a unique 9991 such that 
for t <tQ the following inequality holds 

\'>p2{t,to, ^-20 ,^0) -1p2{t,to, ??2fc_i + V2k-1 (<^01 , /^^O, ato , , ^01 ) I 

+ |'/'2(i,io,^20,¥'o) - ^2{t,to,T^li,_i +V2k-i{^oi,l3to,ato,fi,e),ipoi)\ 
< £3e^''^3(t-to) (1(^0 -(^oil + \^20-^''2k-i-V2k-iimi,Pto,ato,ti,e)\) , (6.12) 

where £3 > 1 and > is independent on a, and e. 

Proof. Let us consider a neighborhood of the point ■02 = '^k- Setting Ci = /3t, (2 = ot, 
new time r = /i^t, and new parameters f?i = /^^, ?/2 = S'^i % = ^"^/l^i X = 1/^^ iii (6-5)- 
(6.6) we change the variables '02 = + V'3 obtain the following autonomous system on 
4-dimensional torus T4 : 

^ = G'{4)^s + G2(03)0i + mSlli^k + ^3, Cl, ^) 

+??2Si2(??^. + 03, ^^2, Ci, C2, /i, e) + ??35i3(??^. + 03, ¥^2, Ci, C2, /i, e), (6.13) 
^ = aoX + 5*2i(t?° + 03,m) + '?i5'2i(i?° + 03, Ci,/^) 

+??2522(t?° + 03, V^2, Cl, C2, Ai, e) + ??3^23(t?° + 03, ^^2 , Cl , C2 , A*, s) , (6.14) 

= PX, -^ = ax, (6.15) 

where G2(V'3)^i := G(??o + V3) - G'(^?0) - G'(^9O)03. 

In the function space C'~^ (T3 x /;^(,) of bounded together with their / — 3 derivatives 
functions w((/32, Ci, C2, A) defined on (i/?2,Ci,C2) G T3, A = (r/i, ??2, m. A*, e) ^ -^Aq = {A : ||A|| < 
Aq}, we introduce the mapping 

Tk{w) = - I exp [-G'{^l)i) Qa{w{(P2(^,C,i^, (2^, A), (p2^, Ci^, (2^, A)^^ 
Jo 

if G'i'&l) > 0, and 

Tk{w) = / exp (-G'(i?^)^) Q4{w{ip2^, C25, A), 9925, C25, A)(i^ 
./ —00 

for G'('!?2) < 0. Here Q4 is the right hand side of (6.13) with exception of the linear term, 
^2ii = V52(C, 9^,(1,(2, A), C15 = /5C + Ci, C25 = a? + C2 is the solution of (6.14) ~ (6.15) for 

03 = W'(¥'2,Cl,C2, A). 

Analogously to the proof of Lemma 5.1 (see also [19]), we apply the fiber contraction 
theorem and show that there exists a unique fixed point 

Wk = ??iVA;i(Ci, A) + mvk2{^2, Ci, C2, A) + rjsvk3{ip2, Ci, C2, A) (6.16) 
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of Tj^{w) in the neighborhood of (0, 0) G C'^^(T3 x Functions in right-hand side of (6.16) 

are C'~^ smooth and 27r-periodic in ip2,C,i,C,2, such that ||t'A:j||c'-* — where positive 

constant M2 does not depend on A, x ^ Xo and a > oq. 

Respectively, there exist /io > and cq > such that for all < ^ < /^o and e < co/i^^^ 

the system (6.13) - (6.15) possesses 2N integral manifolds (6.10). 

The proof of the inequalities (6.11) and (6.12) is analogous to Lemma 6.1 in [19]. ■ 
Corollary 6.2. Under the conditions of lemma 6.1, system (6.1) - (6.2) has 2N integral 

manifolds 

Vj = {{/3t + ^^j + Vj{ip,(3t,at,n,e),ip,t) : ip£Ti,t£R}, (6.17) 
where C''~^ -smooth function 

^3 

Vj = fi'^Vjoi/St, fi) + e'^Vji{(p, l3t,at, fi,e) H Vj2{f,l3t,at,fi,e), (6.18) 

is 2TT-periodic in <f,l3t, and at. Manifolds corresponding to j = 2k, k = 1,...,N, are expo- 
nentially stable for t — t- +cxd and manifolds corresponding to j = 2k — l,k = 1,...,N, are 
exponentially stable for t —00. 

Note that the expressions for the integral manifolds in the nondegenerate case, i.e. gi ^ 0, 
can be simplified. In particular, in this case Sis = and = in (6.5) - (6.6) as well as 
5i3 = and 523 = in system (6.1) - (6.2), since R04 = 0, R14 = R24, = in system (5.6) - 
(5.8). The following corollary gives these expressions. 

Corollary 6.3. Let S13 = S23 = 0. Assume A S [G^,G^] such that (6.9) has 2N nondegen- 
erate solutions ^ = -d^ , j = 1, . . . , 2N. Then there exist fiQ > and cq > such that for all 
< < /^o o,nd < e < Co system (6.1) - (6.2) has 2N integral manifolds 

TZj = {{f3t + ^'^j +Vj{^,f3t,at,n,£),^,t) : (p£Ti,t£R}, (6.19) 

where 

Vj = fi'^Vjo{f3t,n) + e'^Vji{ip2,Pt,at,fi,e), 

with C^~^ smooth, periodic in ip2,(3t,at functions vjk, such that \\vjk\\Qi-2 < M3 with the 
constant M3 independent on a,fi,£. The manifolds 7^2fc; k = 1, ...,N, are exponentially stable 
in the sense of inequality (6.11). The manifolds 7^2fc-i; k = 1, ...,N, are exponentially unstable 
in the sense of inequality (6.12). 

7. Proofs of theorems. We start by proving the degenerate case, since the non-degenerate 
case will be a particular case of this proof. 

7.1. Degenerate case gi = 0. Let us prove Theorem 2.3. In the case gi = 0, two 
averaging transformations (3.1) and (3.7) reduce system (1.1) into (5.1) with = /U^, 

1/a = e^, and g = 52- The obtained system (5.1) is further transformed using the local 
coordinates to (5.6) - (5.8). In the latter system, for any fixed /3 there exists an orbitally 
asymptotically stable integral manifold (5.9) accordingly to Lemma 5.1. Taking into account 
regular dependence of the right-hand side of system (5.6) - (5.8) on /3 G [/3i,/32] we conclude 
that constants fiQ, £q, Mi, k, and L from Lemma 5.1 can be chosen common for all /3 G /32]. 
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Hence, since dependence (3.11), in the original system (1.1) there exists the exponentially 
stable integral manifold 

X = e^^(xo(V') + ^o(^)'") + -MMi^) + ^o(^)^^, at, ^, -), (7.1) 

a a a 

where u is defined by (5.10). Taking into account the properties of functions U2 and u, we 
conclude that this manifold has form (2.17). The equations on the manifold is obtained by 
substitution (5.9) into (5.6) - (5.8) and are given by (6.1) ~ (6.2). Substituting the parameters 
H = 7/\/a and e = l/\/a, one obtains the system on the manifold (2.18) - (2.19). 

All solutions from some neighborhood of the torus Ti are approaching 9Jl2(a, 7). There- 
fore, in order to show the remaining statements of Theorem 2.3, it is enough to consider the 
behavior of solutions on the manifold 9Jt2(a, /3, 7). Note that the manifold 9J{2(o5/9,7) corre- 
sponds to the manifold in the local coordinate system {h, tp, ip) with parameters and 

£. 

We first consider system (6.1) - (6.2) with parameters fi and e. Let us fix any positive e. 
For the set S2 of singular values of G2 we define the following sets 

^2(e) = {5G[G2",G'+]: dist(<7,S2)>e}, 

and 

^2(e) = {0G [0,27r] : ^2(0) G e2(e)}. 

Taking into account that the sets B2{e) and ^2(e) are compact, there exists a positive constant 
<; such that 

dG^ie) ^^^^j^ eeA2{e). (7.2) 



de 

By Corollary 6.2, for fixed A E S2{e), there exist ^0 > and cq > such that for 
all < < //o and e < cq^^^^ system (6.1) - (6.2) has 2N integral manifolds Vj. Since the 
estimate (7.2) is uniform, constants ^0 > and cq > can be chosen common for all A € i3(e), 
and therefore for all a, 7, and /3 satisfying (2.20). Hence, under conditions (2.16) and (2.20), 
there exist 2N integral submanifolds £,j C 9Jt2, j = 1,...,2N of the form (2.21) for system 
(1.1). Locally, the manifolds corresponding to j = 2k, k = 1, . . . ,N are exponentially stable 
and the others are exponentially unstable. 

In order to complete the proof, let us show, that any solution from the manifold dJl2{a, f3, 7) 
is attracted to one of the manifolds £2fc as t — )■ 00. For this, it is more convenient to consider 
the reduced system on the manifold (6.5) - (6.6) and the equivalent to £,j integral manifolds 
Uj (6.10). 

All manifolds Il2kj k = 1, . . . , N , are asymptotically stable with the exponential estimate 
(6.11) and manifolds n2fc_i are asymptotically unstable with the estimate (6.12). Therefore, 
by (6.12), on a finite time interval T, the solution (^2(^)1 9^2 (i)) of (6.5) - (6.6) with initial 
values {'(p2{io),^2{to)) not belonging to the manifold n2fc_i, i.e. 

^2(^0) / i?2fc-i + V2k-i{^2{to), Pto,ato, fi,e), 

and 1^2(^0) — ''^2fc-il ^ '^Oi reaches the boundary of (5o-neighborhood of unperturbed manifold 
'&2k-i ^ '^1' ™oi'e exactly, the value of 'ip2{t) reaches i?2fc-i ~ "^o or '&2k-i ~^ where 60 is 



24 



L. RECKE, A. SAMOILENKO, V. TKACHENKO, S. YANCHUK 



defined from Lemma 6.1. The time interval T depends on the initial values (^^2(^0)) ¥'2(^0)) 
and parameters fi,e. 

Further, any solution {tp2{t) , ^2{t)) starting outside of a 5o-neig hborhood of 1?°^.^ x Ti 
reaches (^o-i^sighborhood of either the manifold '&2k ^ '^1 ^2fc-2 ^ '^1 ^ finite time interval. 
This follows from the fact, that the right hand side of (6.5) is uniformly bounded from zero 
on this set for small enough parameters e and fj, (see also Lemma 6.2, [19]). 

Next, by the inequality (6.11) of Lemma 6.1, any solution from the (5o-neighborhood of 
the manifold i^gfc ^ "^i attracted to the stable integral manifold Il2k- 

As a result, solutions {Tp(t),(p{t)) of system (6.1) - (6.2) that, at initial point t = to do not 
belong to unstable integral manifolds V2k-i, k = 1, N, are attracted for t > to to solutions 
{ijj{t),ip{t)) on one of the stable integral manifolds V2k, i-e. 

tPit) = (3t + + V2k{^{t),lit, at, £), 

and ip{t) satisfies (6.2) with ^(t) = ^(t) and j = 2k 

^ = ao + IJ^S2i{Pt + ^] + Vj,/3t, /i) + eV^5'22(/3t + + vj, 93, /3t, at, ^i, e) 

+ (/3t + 'd^j+Vj,(p, /3t, at, fi, e) . (7.3) 

If a solution {ip{t), (p{t)) of (6.1) - (6.2) at initial point t = to belongs to one of the unstable 
integral manifolds V2k+i then this solution has the following form 

il^it) = /3t + 7?2fc+i + V2k+i{^{t),(3t, at, n, e), 

where ip{t) satisfies (7.3) for j = 2k + 1. 

Lemma 5.1 implies, that any solution {h(t),ip{t),ip{t)) of (5.6) - (5.8) with ||/i(to)|| < 
is attracted to one of the solutions (/i(t), ■0(t), (^(t)) on the integral manifold SCI^^g such that 

h{t) = u{'4){t), (p{t),l3t, at, /i, e), 
^{t) = /3t + ^9° + vj{'f{t),(3t, at, n, e), 
(p{t) is a solution of system (7.3) 

with some j, 1 < j < 2N . Therefore, every solution X2(t) of averaged system (3.10) that at a 
certain moment of time to belong to a small neighborhood of invariant torus Ti tends to one 
of the solutions 

X2(t) = e^^W (xo(/3t + + Vj) + $o(/3t + ^° + Vj)h{t)) , j = 1, . . . , 2iV. 

Taking into account averaging transformation (3.11), and the form (6.18) of Vj, every solution 
x{t) of system (1.1) that at a certain moment of time to belongs to 5-neighborhood of torus 
T2 tends as t — )• +00 to one of the manifolds £j given by (2.21). 

The proof of theorem 2.4 follows from (2.21) taking into account smallness of 7/a, 
and for a, (3, and 7 satisfying (2.20). 
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7.2. Non-degenerate case. If gi ^ 0, the averaging transformation (3.1) in (1.1) leads to 
system (5.1) with 7 = fi^, 1/a = e^, and g = gi, r2 = 0. The existence of the asymptotically 
stable invariant manifold 9Jti(a,/3,7) of the form (2.8) follows directly from Lemma 5.1 and 
the coordinate transformation (5.2). 

For the proof of the other statements of Theorem 2.1 we consider the dynamics on the 
manifold 5!Jti(a, /?, 7). For any 0,7 and /3 satisfying (2.11) there exists an even number of 
points j = 1, 2N, which are solutions of the equation /3 — /3o = jGi{6). The number 
depends on the parameters a, 7 and f3. 

By Corollary 6.3, for all /x G (0, no] and e < cq with sufficiently small /xq and cq system (6.1) 

— (6.2) has 2N integral manifolds TZj corresponding to points Analogously to the proof of 
Theorem 2.3 in the degenerate case, we show that every solution {Tp{t),ip{t)) of system (6.1) 

- (6.2) is attracted for t — )• +00 to some solution {ip{t), (p{t)) on one of the integral manifolds 
TZj, j = 1,...,2N. 

The manifolds TZj in system (6.1) - (6.2) correspond to the manifolds 9Ij in system (1.1). 
Therefore, every solution x{t) of system (1.1) that at a certain moment of time to belong to 
(5-neighborhood of torus T2 tends to one of the manifolds 9Tj as t — t- +00 and satisfies estimate 
(2.13) from Theorem 2.2. 

8. Example 1. In this section, we illustrate the obtained results using a model of a laser 
with saturable absorber with an external electro-optical forcing. The model has the form (1.1) 
with X E M'', 

f fi{a- xi - xi{xl + xl)) \ 



fix) 



^{xi 
^{xi 



X2 - 

- X2 

- X2 



CX2 {xl + xf)) 



.1) 



1) X3, - aoX4 
1) X4 + aoX3 J 

where fi, a, b, c, and oq are real parameters. This system satisfies the symmetry condition 
(1.2) with matrix 








A 



/o 




Vo 



-1 



^.2) 



/ 



Note, that in the unperturbed case, the model can be written in partially complex form 



'I I |2 

Xi = fj, [a — xi — xi \y\ 



x'2 = fJ. (b - X2 - CX2 |yn , 



y' = ^{xi - X2- l)y + iaoy, 



(8.3) 
(8.4) 
(8.5) 



where y = X3 + ix^ , i = \/— T, which can be further reduced to the Yamada model describing 
the dynamics of the intensity of a laser with saturable absorber [27, 7]. In this case, y 
plays the role of the electric field with intensity / = |yp = x"^ + x\, and x\, X2 are gain 
and absorption, respectively. In [7], the parameter regions are numerically obtained, for 
which the Yamada model (in coordinates {xi,X2,I)) has an asymptotically stable periodic 
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solution (^Xi{f3ot) , X2{(3ot) , {f3Qt)) , I^{/3ot) > 0. We choose the parameter values a = 7.0, 
b = 5.8, c = 1.8, and /x = 0.04 belonging to this region. For these parameters, /3o = 0.071 
and the corresponding complexified system (8.3) - (8.5) has the quasiperiodic solution x{t) = 

^a;5(,5ot), X2(/3ot), y^/^^otje*^""*^^^^ , where is an arbitrary real constant resulting from the 

symmetry. Note that the S^- equivariance of the right-hand side of (8.3) - (8.5) with respect 
to the transformation y — t- e^^y is equivalent to the S'^-equivariance (1.2) of the function (8.1) 
with matrix A from (8.2). The corresponding solution (1.5) of (1.4) with / defined by (8.1) 
reads then 

x{t) = (x;(/3ot), xO(/3ot), V/°(/3ot)cos (aot + , V^'^Wot) sin {aot + O)^ ■ 

In this way, using numerical results from [7], we see that the unperturbed system (1.4) 
with (8.1) has an exponentially attracting invariant torus (1.7) with A from (8.2) and 

xoW = [x'ii^P),xU^),xl{^P) = ^^1^,0)^. (8.6) 

Here xq{^) is the periodic solution of the system (2.1) transformed to the corotating coordi- 
nates. It can be found numerically using direct integration. Figure 8.1 illustrates xq{iIj) for 
the chosen parameter values. It has shape of a pulse typical for lasers with saturable absorber. 

Accordingly to Theorems 2.2 and 2.4, the synchronization region for the parameters a, /3, 
and 7 is determined by the properties of function Gi{4>) from (2.6) in the nondegenerate case 
and G2(V') from (2.15) in the degenerate case. In both cases, one needs the periodic solution 
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Figure 8.2. Periodic solution p\ (tp) of the adjoint variational system (2.4) for unperturbed Yamada model. 
The eomponents pi,i, Pi,2, and pi,3 are shown, while pi^i^ip) = 0. 



pi {tp) of the adjoint variational equation (2.4), which satisfies the orthonormality condition 
(2.5). Using standard numerical methods, which involve computation of the monodromy 
matrix, we found pi (ip) numerically. Because of (2.5), pi is defined uniquely. The result is 
shown in Figure 8.2. 

With the given pi (tp), the effects of arbitrary perturbation of the form (1.1) ~ (1-3) can 
be studied. For example, let us consider the perturbation 

75 (x, I3t, at) = 7 {gel (Pt) , 0, cos (at) gop {Pi) , sin {at) gop (/3t))'^ . (8.7) 

This kind of perturbation may correspond to some electro-optical external injection, where 
^gd {at) stands for an electric and 7 cos {at) gop {j3t), 7 sin {at) gop {(it) for optical components, 
respectively. Here, we consider (8.7) just as an illustrative example for our method. 
The function gi{x,ip) from (1.9) is then reduced to 

gi (x,V) = bci(^), 0,0,0]^ 

and ^ 

Gi {'4^) = ^ jj gox{e-ij)pi,i{e)de. 

Accordingly to Theorems 2.1 and 2.2, the locking regions have the form as in Fig. 2.1 and the 
corresponding angles are determined by the extrema of Gi (V'). Let us consider, for example, 
Oci (V') = sin ('(/'). Then the function Gi{ip) has the form 

Gi {ip) = Gs sin ij: + Gc cos ip 
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with 

I /'27r -1 rln 

Gc = -— pi,i{9)cosede, Gs = — pi,i{e)smede. 

^71" Jo Jo 

For the chosen parameter values (using pi{ip) given in Fig. 8.2), we have Gc ~ 1.67 and 
Gs ~ —2.33. Hence, Gi (ip) is a harmonic function with extrema Gf = sj G^^ G\ ~ 2.86 
and G^ = —G\^ see Fig. 8.3(a). The synchronization domain is symmetric and given by the 
conditions a > qi, < 7 < 71, and |/3 — /3o| < 7^^ w 2.867. The bounds 71 and ai cannot 
be determined using our approach. If parameters belongs to synchronization domain, then 
for any solution x{t) of system (1.1) with right-hand side defined by (8.1) and (8.7) which is 
at a certain moment close to invariant torus (1.7) of unperturbed system with xq from (8.6), 
there exists a € M such that 

\xi{t) - x?(/3i + a)\ + \x2{t) - x^i^t + a)\ + \xl{t) + xl{t) - I^{^t + a)\ ^ 

for large t. 

Interestingly, when the electric perturbation has two maxima, e.g. gei{ip) = 0.5 cos + 
sm(2ijj), then Gi {ip) can also have two maxima (see Fig. 8.3(b)) which leads to the appearance 
of multiple stable synchronized manifolds (Theorem 2.1). In this case, the set of singular 
values is 5i = {G^ = 1.87,1.22,0.01,-2.77 = Gj^} and the synchronization domain is given 
by (2.11). 

In a similar way, the influence of other perturbations of the form (1.1) - (1-3) can be 
studied. 
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9. Example 2. In the next example, the system has quasiperiodic solution of the form 
(1.5), which can be written in simple analytical form. We consider system (1.1) with x = 

(xi, X2, xs, X4) G and 



(1 



I3qX2 + xi{x1 + xl- 

-/3oXi + X2ixl + X4 



XT 



X^)X3 + (/3oXl + /3oX2 + Qo)x4 



\ -{l3oXi + /3oX2 + ao)x3 + (1 



xf 



xi)x4 j 



g{x,ip,ip) 



/ X3 cos (f + X4 sin 93 \ 
—2:3 sin ip + X4 cos ip 
(X3 + xl) sin ip cos 99 
\ (xl + xl) sint/'sinc^ / 



This system satisfies symmetry conditions (1.2) and (1.3) with matrix A defined by 
Unperturbed system (if 7 = 0) has the solution 



.2). 



x{t) 



o*xo(/3ot) 



\ 



sin(/3oi) 
cos(/3ot) 
sin (sin(/3ot) — cos(/3ot)) 
\ cos (sin(/3ot) — cos(/3oi)) / 



xO(/3ot) 
V 2:4(^0^) / 



The corresponding variational equation (2.2) has two linear independent periodic solution 

dxo(V') 



91 (V') 



dtp 



92 (-0) = -4x0 (V')- 



The adjoint equation (2.4) has two periodic solution pi{ip) and P2{ip) such that pT('i/')'?fc('0) 



6jk for all ^ and j,k = 1,2. It can be verified that Pi{iIj) = {cosip 
computation we find that gi = and the second averaging function 



sin-;/', 0, 0). By direct 

fl'2(a;o(V'),'0) = (0, -sin-i/', -xl{tp)sm'^ ip,xl{'4))sm^ tp) 



and 



1 /" 

G2W = ^ J Plii^ + 0)92ixoii^ + 0),e)de 
1 /■^'^ 1 

= — / sm{Tp + 6) sin = — cos ip. 
2vr Jo 2 



Therefore, Gg" 



1/2 and by Theorem 2.3 the synchronization domain (where 



condition (2.22) is fulfilled) is defined by |/3 - /3o| < j'^G^/a 
ci/q < 7 < C2\/a with some positive constants ai, ci, and C2. 



J /2a for a > ai and 
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